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Abstract 

Searching for a parking spot can waste time and gasoline. This waste can be reduced by 
assigning drivers to parking lots based on their destination and arrival time. In such a system, 
drivers could request a parking spot in advance and be alerted (e.g., via their phone or vehicle) of 
their assignment to a specific parking lot or available spot. In this paper, a parking assignment 
system is described to allocate parking spaces in a fair and equitable manner. Heuristics are 
developed to solve the underlying large scale optimization problem. The efficacy of the system 
is demonstrated by applying our algorithms to real data sets. 


1 Introduction 

A huge societal problem arises from drivers searching for a parking space. Statistics illustrating the 
waste associated with this issue border on the unbelievable. For example, in the United Kingdom, 
the average person spends 2,549 hours looking for a parking spot over their lifetime [T]. Further, a 
survey conducted by Allianz Insurance estimated that 95% of British people added an extra mile 
to their trip searching for a parking spot [2]. Such statistics arise elsewhere and are not just in 
the UK. For example, it was recently reported that over one year in a small Los Angeles business 
district, cars cruising for parking burned 47,000 gallons of gasoline and produced 730 tons of carbon 
dioxide [3]. Meanwhile, the consulting firm McKinsey recently claimed that the average car owner 
in Paris spends four years of his or her life searching for a parking space [4] . The parking assignment 
problem associated with electric vehicles (EVs) becomes even more acute. Due to the limited range 
of these vehicles, the marginal cost of expending energy to search for spaces may, in some cities, 
be prohibitively high. Thus, there is a real and compelling societal and economic need to revisit 
parking. 

Currently, parking guidance systems can alert drivers to the availability of parking spots in 
given areas [5]. However, this type of system may not work in a crowded area with few available 
spots. Moreover, if parking spaces are not specifically assigned to vehicles, then a spot may no 
longer be available by the time the vehicle arrives. Assigning drivers to parking lots could reduce 
the time and gasoline wasted by drivers searching for parking spots and could ensure that the 
allocation of drivers to parking spots meets specified criteria. 

An assignment of drivers to parking lots can be seen as a resource allocation problem, where 
the resource to be distributed is parking spots. An important quality of resource allocation is 

*IBM Research Ireland, Damastown Industrial Estate, Dublin 15, Ireland. Email: nicole.taheri@ie.ibm.com. 

The authors would like to thank TAPAS Cologne for letting us use their data. 


I 



fairness; the assignment of goods to users should be fair given a specified definition of fairness [6]. 
For a parking lot assignment, fairness can be measured in terms of the travel time between user 
destinations and assigned parking lots. 

Following on from our prior work in this area [am, we propose an optimization-based method 
to assign vehicles to parking spots in a fair way. The principal conceptual difference between our 
previous work and the work presented here is that we consider a centralized system that solves 
an exact optimization problem (with some relaxation), and is built on an elaborate notion of 
fairness from the viewpoint of the vehicle owner. Specifically, we define a driver-centric measure 
to determine the fairness of a parking lot allocation, and construct an algorithm to optimize this 
fairness measure. We then test our method on two sets of real driving data: one from Cologne, 
Germany mEU] and one from New York City taxis mm- Our results show that using our 
method could improve the fairness of a parking assignment by around 140% when compared to 
basic methods with no optimization. 

1.1 Related Work 

Within the research community, questions concerning how to manage parking space supply-demand 
mismatch are actively being investigated from a variety of different angles. One important aspect 
concerns the delivery of up-to-date, accurate, real-time information to parking systems in order 
to achieve greatest system efficiency. Parking guidance and information (PCI) systems can help 
inform drivers of where there are available spots; for example, a system could list the number 
of available spots on each floor of a parking garage or in each parking lot in a city [5]. Such a 
system requires detecting vehicles in parking spots, and there have been a number of works that 
consider the technology (i.e., hardware and/or software) required to detect vehicles in parking 
spaces, including HMS]. The effects of parking guidance on a city have also been analyzed mM- 
However, these systems do not specifically assign vehicles to parking spots, they only make drivers 
aware of the availability in the specified areas. This type of system is not ideal for street parking, 
or situations when there are only a few parking spaces available, as the availability of a parking 
space in a crowded area can change quickly. This type of system could also result in unnecessary 
overcrowding of certain areas 

A predictive approach is proposed by Caliskan et al. m and further studied in [20] . The authors 
develop a method to predict the likelihood of a parking space being available at the estimated 
time that the car will arrive there. In m, Ceng et al. view the parking problem as a dynamic 
resource-allocation problem. In this paper, algorithms are developed to assign drivers to parking 
lots dynamically as they leave their destinations. The objective of their formulation is to minimize 
the sum of the utility functions for all the drivers. Chou et al. [22| use a different model and 
consider a network of parking lots with negotiable prices; parking lots are selected for drivers in 
a way that benefits both the drivers and the car park operators. The work of Teodorovic and 
Lucic |23] proposes a method to define rules for assigning drivers to parking spots using fuzzy logic 
and integer programming techniques. The parking assignment problem has also been solved using 
a stable matching algorithm [23], which yields a Nash equilibrium solution, and through differential 
pricing |25j . Parking lots can also be assigned to drivers through matching under preferences [26]. 
Griggs et al. created a distributed privacy-preserving scheme to allocate parking efficiently with 
quality of service guarantees [8], where fairness is one such objective. Finally, in [7] cars are assigned 
to lots in a balanced manner, where fairness is considered from the view point of the parking lot 
owner. Our work follows the prior work described in [7118] . However, there are important differences. 
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We provide an advanced driver-centric notion of fairness, and give a heuristic that solves a very 
complicated optimization problem that goes far beyond that described in any of the aforementioned 
references. Importantly, our fairness notion can be extended to include other aspects of fairness 
(aspects associated with EV ownership) as a particular need arises. 

2 Problem Statement 

We propose an optimization-based method to assign parking lots to drivers in a fair way. We 
assume that the starting and ending times of all trips are known at the beginning of the day (i.e., 
a 24-hour period), as well as the starting and ending locations. 

Without loss of generality, we assume that a vehicle can be collected from a parking lot as well 
as dropped off (i.e., a driver may pick up their vehicle from a specified parking lot and park it at a 
different parking lot near his/her destination). We assume that a driver will pick up a vehicle from 
the parking lot closest to their origin location and our method will assign the driver to a parking 
lot in which to drop off the vehicle that is near their destination. 

2.1 Notation 

Table [T] defines the terminology to be used in the remainder of the work. 

Table 1: Sets, Constants and Variables 


Sets 


Name 

Description 

7^ 

Drivers 

r 

Time Periods 

£ 

Parking Lots 


Constants 

Name 

Description 

“M € 1R2x|£| 

X € 

s,d elR2x|7e| 

cx G 

Location of Parking Lots 

Maximum Number of Parking Spots in each Lot 

Start and End Location of Driver Trips 

Start and End Time of Driver Trips 

Average Walking Speed (in miles per hour) 


Variables 


Name 

Description 

y 


Destination Lot Assignment 

x{t) 

€ 


Number of Filled Parking Spots at time t 

Y{t) 

€ 


Number of drivers arriving at each lot at time t 

Z{t) 



Number of drivers leaving each lot at time t 

D 

e laf' 


Distance from Assigned Parking Lot to Destination (in miles) 

/3 



Total Overhead Time for each Driver Trip (in hours) 


3 













For all trips r G 7^ and parking lots £ € C, we define the binary decision variable y G {0,1}^^^ 
and the binary parameter z G {0,The parameter zir will specify from which parking lot 
the vehicle will depart (its origin), and the variable yir will specify in which parking lot the driver 
will park (its destination). Because we assume that each vehicle will start in the parking lot closest 
to the driver’s origin location, z is a parameter, not a variable. However, the lot in which a vehicle 
will park near the destination is a decision variable to be optimized. The values of y and 2 ; can be 
defined with the equations: 


Vlr — 



if trip r ends up in lot £ 
otherwise 



if trip r starts from lot £ 
otherwise. 


The values of 2 ; can easily be calculated given the origin locations of the trips, s. 

We also define the state variables Y,Z^ gl^l^lTI that keep track of the number of vehicles 
arriving and leaving each parking lot over time: 

Yl(t) — ^ ^ yir) ^ ^ Zlri 

r£7l:t^=t re'R:t^=t 


where are the origin and destination locations of driver r respectively. The variable xi{t) 

represents the number of filled (or parked-in) spots in lot £ at time t, and is the maximum 
number of available spots in parking lot £. A feasible set of lot assignments and parking lot 
states xi{t) will be in the set: 


n = 


< {x,y) : 


xe{t) = Xi{t - 1) + Yi{t) - Zi{t) 
xeit) < xi 

lec 

Vir G {0, 1} 


w G £,t G r' 
W eC,t£T 

Wen ' 

WeC,ren 

/ 


(FEAS) 


The first equation in (IFEASI) updates the state of each parking lot over time: the number of 
occupied parking spots is equal to the number of occupied spots in the previous time period, plus 
the arriving vehicles and minus the departing vehicles in the previous period. The second equation 
in (jEEASI) ensures the number of vehicles in each lot is less than the capacity, the third equation 
makes sure that each trip ends in only one parking lot, and the last equation constrains y to be 
binary. 


3 Optimizing Fairness 

Assigning drivers to parking lots leads us to the following questions: Is there a fair way to assign 
drivers to parking lots? And what does it mean for an assignment of parking lots to be fair? 

In our problem, a fair assignment will be one in which the distance between drivers’ destinations 
and their assigned parking lots are similar. Let fir be the amount of time driver r e n spends 
traveling from the assigned parking lot to his/her destination: 


fir = Oi 


dr 


'^{yir ■ Mi) 
£ec 


( 1 ) 
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We define fairness as a lack of envy among drivers, where the envy between two drivers ri,r 2 G 7^ 
is defined as: 

Eri,r2 l/^ri “ /3r2 I ) (2) 

where j3r is the walking time for driver r ^ TZ from his/her parking lot to destination. 

We define fairness as the mean value of £lri,r 2 over all pairs of drivers ri,r 2 G TZ, where a 
smaller value is better. To attain an assignment of drivers to parking lots that optimizes this 
fairness measure, we can minimize the mean of the envy values over all pairs of drivers with the 
objective function: 


F{f5) = mean (£1^1,r 2 ) = (1/ l-R/ 

ri,r2&'R- 


ri,r2€'R. 


\f3ri - A 


T2 I 


(3) 


Definition 1 A fair assignment of destination lots to drivers results in a value of F{I3) that is as 
small as possible. 

3.1 Formulation 

We assume that we are given the origin and destination locations and times for each driver, 
{s,d,F,t^) at the start of the day (i.e., a 24-hour period). Let g : ^ 1R be a function 

that takes the waiting time of each driver r G 7^ and assesses the fairness, where a smaller value 
implies a more fair assignment. Then we can formulate the problem that maximizes fairness as: 


minimize g{fd) 
subject to (x, y) G n 


jdr = a 


dr 


• Me) 

iec 1 


(OPT-FAIR) 


Vr G 7e (4) 


The first constraint ensures that the lot assignments defined by x and y are feasible for the 
given drivers and parking lots. Constraint (HD calculates the distance from the driver’s destination 
to the assigned parking lot and multiplies it by the average walking speed, a, to get the total travel 
time for driver r. The formulation (lOPT-FAIBll is a mixed-integer linear program (MILP) that 
can be solved in CPLEX |27j . 

This optimization problem can be solved for a chosen time period (e.g, a single day); our aim 
is to choose an objective function g{f3) for (lOPT-FAIRI) . so that the average overhead waiting time 
over a set of drivers is fair. 


3.2 Intractability of F(/5) as the objective 

Using the definition of fairness based on the function F(/3) in ([3D, an optimally fair lot assignment 
can be found by setting the objective function in (jOPT-FAIBI) to g(/3) = F{j3). That is, ideally 
we could minimize the mean of all the envy values, or differences between driver walking times. 
However, setting y(A) = F{j3) as the objective in (lOPT-FAIRp results in an intractable optimization 
problem; taking the absolute value of each pairwise difference in travel times results in a large 
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problem that cannot be solved in reasonable time. For example, we used CPLEX to try to solve 
the problem (lOPT-FAIEI) with the objective g{l3) = F(/3), for |£| = 10 parking lots, \TV\ = 100 
drivers, and |T| = 24 time periods, and the solver did not even find a feasible solution after 12 
hours. 


3.3 Choosing an Objective Function 

Given the intractability of using F(/3) as the objective, we would like to construct a method that 
achieves a similar result, but that is computationally tractable. The idea of minimizing F(/3) in 
(|OPT-FAIPp is to obtain a lot assignment in which all drivers’ walking times are similar. We 
propose a method to find a solution of (jOPT-FAIRp that achieves this same goal, but with an 
objective that makes this problem tractable. 

Given a feasible set of lot assignments y € {0, Ijl^l^l^l with corresponding walking times /3 G 
the mean walking time is defined as: 

H{P)-.= {1/\TZ\)Y,^. (5) 

rdTZ 


A fair assignment of drivers to parking lots can be attained by minimizing the absolute value of the 
differences between each driver’s walking time and the mean walking time, H(/3). This will provide 
a solution where the walking times of the drivers ares similar and as close to the mean value as 
possible. 

Given a feasible solution (x, y, /3) to (jOPT-FAIRp . the mean walking time H (/3) can be calculated 
and the problem (lOPT-FAIPp can be re-solved with the objective: 


G{f^) = ^\Pr-H0) 

rGlZ 


( 6 ) 


Our method is based on solving (lOPT-FAIRI) with this objective, which will yield a tractable 
problem and a fair solution. 

To further reduce the size of the problem, we can consider only a subset of the lot assignments 
as decision variables. That is, because some of the walking times Pr may be sufficiently close to 
the mean value H(f3) already, we can fix the lot assignments within a chosen range of H{I3) and 
only change those with a large or small travel time. Define the set S ^TZ: 

S := {r G 7^ : € [(1 - e)H0), (1 + e)/7(^)]}, (7) 


where e > 0 is a chosen value that determines the desired range of walking times from the mean. 

For a feasible set of lot assignments, y G {0, with corresponding walking times (3 G IRI^I, 

we define the set of constraints based on S as: 




{x,y) G D 

' 

h{s) = < 

{x,y,P) 

/3r = Q: ■ II dr ~ ' ^i) ||]^ 

Uir — Vir 

Vr G 7^ 
yreS,ie£ 


K 

Pr = A 

Vr G 5 


and we define the new objective function: 

G{P,S) = '£\l3r-H0) 

r^S 
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The following optimization problem will then find new lot assignments that optimizes fairness and 
reduces the envies among drivers’ walking times: 

minimize G(/3, S) (OPT-MEAN) 

subject to {x,y,P) € Q{S) 

Notice that after sol vine the MILP (jOPT-MEANp and obtainine a new assienment (v*. 
the mean walking time will change according to the new lot assignments and walking times /3*. That 
is, the mean walking of the new solution, i7(/3*), will differ from H(P), and thus (jOPT-MEANjl 
can be re-solved with the new mean walking time, H(/3*). This process of re-solving can continue 
until the mean walking value converges (i.e., stops changing). Iteratively re-solving (lOPT-MEANp 
with the updated mean value H{j3), and the corresponding updated set S will reduce the envies of 
all the drivers until their values cannot be any closer. 

Our full method is described in Algorithm [TJ This method iteratively minimizes the difference 
between a subset of the drivers’ walking times and the mean walking time of the previous iteration, 
until the solution has converged. The algorithm takes three parameters: e is the desired range 
from the mean walking time, 6 is the tolerance used to determine if the algorithm has converged, 
and maxiter is the maximum number of iterations. 

Algorithm 1 Minimizing the envy 
1: Choose values for €,S, maxiter 

2; Find a feasible lot assignment € {0,with corresponding walking times G 
3: i <— 1, converged ^ 0 
4; while not converged do 

5: S^^{re7Z: € [(1 - (1 + e)i/(/3*-^)]} 

6: (x\y^,/3^)<— argmin G(/3,5*) 

(a;,y,/3)60(5b 

7: if lff(/3^) — H{13'' ^)| < (5 or z > maxiter then 

8: converged ^ 1 

9: end if 

10: i ^ (i + l) 

11: end while 


4 Results 

We tested Algorithm [T] on two sets of real driving data. A description of the data and the results 
are below. We also compared our method to the parking assignment formulation considered by 
Geng and Cassandras in [21[I28| . All of our implementation was done in Matlab [29] and used IBM 
ILOG CPLEX |27j to solve the optimization problems. 

4.1 Data 

Real trip data was collected from New York City taxis in 2013 imiiii. From this dataset, we 
extracted information on the origin and destination locations and times for over 50,000 trips in 
January 2013. 
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We also used real driving data from Cologne, Germany, collected for the Simulation of Urban 
Mobility (SUMO) project pfTn]. This data is known as the TAPAS Cologne data and provides us 
with driver behaviors in Cologne for a 2-hour period. 

In order to model a system with known parking lot locations, for each data set we created 10 
parking lots in the given area. For each set of drivers and their known destinations d € we 

used A:-means clustering m on the locations d to hnd 10 clusters of destinations; we then set the 
locations of the 10 parking lots to the centroids of each of the clusters. 

Let w ~ U[a,b] signify that u; is a pseudorandom integer from a discrete uniform distribution 
between o and b. For each of the parking lots, we assigned a capacity and initial occupancy 
Xi{0) of: 


X, um\/\c\) + i, {\n\/\c\) + 2] 

Xi{0) ~ W[(x£/4), (3-X£/4)]. 

In other words, the capacity of each parking lot will be the average number of parking spots needed 
in each lot (|7^|/|T|) plus either 1 or 2, and each lot will initially be between one-quarter and 
three-quarters full. The average walking speed a was assumed to be 5 km/hour |31] . 

4.2 Fairness Results 

To determine the relative fairness of our method, we also implemented the two methods described in 
the algorithms below. Algorithm [2] (Minimum Sum) solves the optimization problem ()OPT-FAIE| 
once with the objective of minimizing the sum of the travel times between the parking lot and 
destination. Algorithm [3] (No Scheme) simulates what might happen with no scheme, if drivers 
were to park in the closest available lot to their destination. (Note that this algorithm will perform 
slightly better than no scheme, because it assumes drivers know exactly where the best available 
spot is.) 


Algorithm 2 Minimizing the sum 
1: Find the lot assignments by solving the optimization problem: 


{x, y, /3) ^argmin ^ I3r 

{x,y)&n 


subject to jdr = ot- 


dj' 


■ Mi) 

i&C 1 


Vr e 


Algorithm 3 No scheme 
1: for t = 1,..., |r| do 

2: for all drivers r € 7^ arriving at their destination at time t do 

3: Find the closest lot to the destination of driver r that has spaces available and assign the 

driver to this parking lot 

4: end for 

5: end for 
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The figures below show a comparison of the fairness of our algorithm to the fairness of Algo¬ 
rithms [2] and [3l For the data for New York City, we ran 500 tests on the data, where each test 
randomly extracted a different set of 100 drivers over a 24-hour period from the data set described 
in Section Ol For each run, we found the mean envy value for that run (i.e., we calculated the 
value of F[j5) in ([3|) for the solution to each method). 





(a) Progression of F{j3) over iterations of Algorithm [T] (b) Probability of F(/3) values 


Figure 1: Implementation Results on NYC Taxi Data 


Figure [Tal shows a progression of the mean envy values, F(/3), over iterations of Algorithm [T] 
in comparison to Algorithm [2] (Minimum Sum) and Algorithm [3] (No Scheme). Figure flbl shows 
the probability that the mean envy F(/3) will be greater than a chosen value. That is, the x-axis 
in Figure [Tb] corresponds to a given mean envy value 7 , and the y-axis is the probability that the 
mean envy F{j3) is greater than 7 for each method. For example. Figure [Tb] shows that 60% of 
envy values for Algorithm [3|(No Scheme) are greater than 4.79 minutes, and 60% of envy values for 
Algorithm [1] (Minimum Envy) are greater than 3.89 minutes. Moreover, the mean improvement of 
Algorithm [ 1 ] over Algorithm [2| is 28.4%, and the mean improvement of Algorithm [ 1 ] over Algorithm 
Elis 25.6%. 

The same tests were run on the data from Cologne, described in Section 14.31 We ran 500 tests 
on the data, where each test used randomly extracted driving data for 100 drivers over a 2 -hour 
period and the data was extrapolated to cover a 24-hour period. 

Figures [2a| and [2b] show the results for the Cologne data that are analogous to Figures [la] and 
llbl for the New York City data. Figure I2bl shows that 60% of envy values for Algorithm [3] (No 
Scheme) are greater than 2.89 minutes, and 60% of envy values for Algorithm [1] (Minimum Envy) 
are greater than 1.23 minutes. Moreover, the mean improvement of Algorithm [T] over Algorithm [ 2 ] 
is 139.6%, and the mean improvement of Algorithm [ 1 ] over Algorithm [3] is 100 . 2 %. 
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(a) Progression of F{j3) over iterations of Algorithm^ 





(b) Probability of F{li) values 


Figure 2: Implementation Results on Cologne Driving Data 


The Jain’s fairness measure is commonly used in congestion control protocols to measure fairness 
of a resource allocation [6ll32ll33j . The value of the Jain’s measure is between 0 and 1, where a more 
fair assignment results in a higher value. For a set of travel times Pr, the Jain’s fairness measure is 
defined as: 


J{P) 


I’^l ■ J2r£7l f^r 


(9) 


We computed the Jain’s fairness measure for each dataset for each method; Figure [3] shows the 
probability that the Jain’s measure J{l3) will be greater than a chosen value. 

Note that, as opposed to in Figures ()2b(l and (llbjl . for the Jain’s fairness measure, it is better 
for the values to be higher. 

For Figure [3] the x-axis corresponds to a given value 7 , and the y-axis is the probability that 
the Jain’s measure c7(/3) is greater than 7 for each method. For example, using the Cologne data. 
Figure [Tbl shows that 60% of Jain’s measure values for Algorithm [3] (No Scheme) are greater than 
0.89, and 60% of Jain’s measure values for Algorithm [1] (Minimum Envy) are greater than 0.97. 
Using the data from NYC taxi trips, 60% of Jain’s measure values for Algorithm [3] (No Scheme) are 
greater than 0.86, and 60% of Jain’s measure values for Algorithm [T] (Minimum Envy) are greater 
than 0.91. 

Moreover, using the Cologne data, the mean improvement of Algorithm [T] over Algorithm [2] is 
8.95%, and the mean improvement of Algorithm [ 1 ] over Algorithm [3] is 7.67%. Using the NYC taxi 
trips, the mean improvement of Algorithm [T] over Algorithm [ 2 ] is 6 %, and the mean improvement 
of Algorithm [ 1 ] over Algorithm [3] is 4.61%. 
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(a) Jain’s measures using NYC Taxi Data (b) Jain’s measures using Cologne Driving Data 


Figure 3: Probability of Jain’s Measure, c7(/5) 


4.3 Comparison to Other Methods 

We compared our algorithm to the method used by Geng and Cassandras in [2T1[28] . In this work, 
the authors assign parking lots to drivers dynamically as they leave their destinations (that is, 
the driving patterns are not known ahead of time). Geng et al. do not consider fairness of lot 
assignments; their method minimizes the sum of the utility functions of each driver requesting a 
parking spot at a given time. 

We implemented the method of Geng et al. as described in the paper [2T] and compared their 
method to ours. We generated random data using the same methods they described, according to 
the specified distributions and parameters meant to simulate parking around Boston University. 
We also implemented our method on a dynamic basis, so that our algorithm was implemented to 
assign parking lots to the drivers that were currently requesting parking spots. 

Geng et al. dehne the utility function of each driver r C 7^ as 


T — \ _L n _ X 

~ Mr + ^ Dr '■ 


( 10 ) 


where Xr G [0,1] is a weight defined by the user, Mir is the monetary cost for driver r to park in 
parking lot £, and D^r is the distance from driver r’s destination to parking lot i. The values Mr 
and Dr are also provided by the user, and these are the maximum values of M^r and D^r that user 
r is willing to accept. 

For each time step k = 1,2,... , let P{k) C 7^ be the set of drivers that have requested a 
parking spot after time k and have not yet arrived at their destination by time k. Using notation 
consistent with Table [T] and adding the time step parameter k, Geng et al. find a lot assignment 
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y{k) G {0,by solving the problem 


minimize y£r{k) ■ J^r + \^~'y^yt.r{k)\ (SMART-PARK) 

£^C,r&V(k) rev{k) V l£C / 


subject to X£{t) = X£{t — 1) -|- Y£{t) — Z£(t) 

yieC,t€T 


X£{t) < X£ 

yieC,t£T 


'^yir{k) < 1 

Vr G V{k) 

(11) 

££C 



y£r G {0,1} 

W€C,reV{k) 


^ ^,{y£r{k^ • J£r) Y 'y ^jiy£r{k 1) ' J£r) 


(12) 


l&C l&C 


Notice that the only difference between the hrst four constraints in (jSMART-PAEK]) and the set de¬ 
fined by n in (IFEASp is the equation dill) , where this equation is an inequality in dSMART-PARK]l 
and an equality in dFEASp . This implies that each driver may not necessarily be assigned to a 
parking spot in dSMAPT-PARKl) , but is guaranteed a parking spot assignment with the constraints 
in (]FEASp . The last constraint dl2p ensures that any new assignment to driver r is no worse than 
their assignment from the previous period, k — 1. 

To compare our method, we used the formulation dSMART-PARK]) . but changed the objective 
function to: 

minimize V {y£r{k) ■ Jtr - J{k)) + y] 1 -'^Virik) , (13) 

^ l&C,r&V{k) r£V{k) V / 

where 

Yp(W\ { X/ ytr{k - 1) ■ Jirj ■ 

' ^ \£ec,rev[k) ) 

This is analogous to the objective defined in ([6]), where the only difference is that (fT^ substitutes 
the travel times /3 with the utility function values J. We then use the iterative method defined in 
Algorithm [T] with the optimization problem defined by the objective (jl3p and the constraints to 
dSMART-PARK]) . This allows us to compare the fairness of our algorithm with that of Geng et al. 

Table [2] compares the relative improvement in the mean of the envy values E(/3) as defined in 
(|3p, over 500 runs. Each run uses different randomly generated data using the guidelines provided 
in m, and compares the fairness of Geng et al.’s method to Algorithm [T] using the formulation 
(I SMART-PARK I) with the objective (fT^ . The values in Table [2] are the mean of the relative 
improvement in the fairness measures E(/3) and i7(/3) for the solutions of each of the 500 runs. 


Mean Envy 

Jain’s Eairness Measure 

14.51% 

4.56% 


Table 2: Mean Relative Improvement in Pair ness using Algorithm [T] over (ISMART-PARKll 
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5 Conclusion 


In this work, we constructed a measure to determine the fairness of a given assignment of drivers to 
parking lots. We used this measure to construct an optimization-based algorithm that iteratively 
reduces the envy among drivers (or iteratively increases fairness) and finds a fair parking lot allo¬ 
cation. Our method was tested on real driving data and compared to a number of similar methods 
to show that our algorithm provides fair parking lot assignments. 
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